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Radar has been proposed as a way to track wake vortices to reduce aircraft spacing and 
tests have revealed radar echoes from aircraft wakes in clear air. The results are always 
interpreted qualitatively using Tatarski’s theory of weak scattering by isotropic atmo- 
spheric turbulence. The goal of the present work was to predict the value of the radar 
cross-section (RCS) using simpler models. This is accomplished in two steps. (1) First, 
the refractive index is obtained. Since the structure of aircraft wakes is different from at- 
mospheric turbulence, three simple mechanisms specific to vortex wakes are considered: 
(A) Radial density gradient in a 2D vortex (B) 3D fluctuations in the vortex cores* (C) 
Adiabatic transport of atmospheric fluid in a 2D oval surrounding the pair of vortices. 
The index of refraction is obtained more precisely for the 2D mechanisms than for the 
3D ones. (2) In the second step, knowing the index of refraction, a scattering analysis is 
performed. Tatarski’s weak scattering approximation is kept but the usual assumptions 
of a far-fieid and a uniform incident wave are dropped. Neither assumption is generally 
valid for a wake that is coherent across the radar beam. For analytical insight, a simpler 
approximation that invokes, in addition to weak scattering, the far-field and wide cylin- 
drical beam assumptions, is also developed and compared with the more general analysis. 
The predicted RCS values for the oval surrounding the vortices (mechanism C) agree with 
the experiments of Gilson (1992) conducted over a wide range of frequencies. However, 
the predictions have a cut-off away from normal incidence which is not present in the 
measurements. Estimates suggest that this is due to turbulence in the baroclinic vorticity 
generated at the boundary of the oval. The reflectivity of a vortex itself (mechanism A) 
is comparable to that of the oval (mechanism C) but cuts-off at frequencies lower than 
those considered in all the experiments to date. The RCS of a vortex happens to peak at 
the frequency (about 49 MHz) where atmospheric radars (known as ST radars) operate 
and so the present prediction could be verified in the future. Finally, we suggest that hot 
engine exhaust could increase RCS by 40 db and reveal vortex circulation, provided its 
mixing with the surroundings is prevented in the laminarising flow of the vortices. 


1- Introduction 


1.1. Motivation 

The photographs of fluid motion in Van Dyke’s (1982) album are alt made possible by 
an optical effect in the visual range of the electromagnetic spectrum, for instance light 
scattering by smoke and the schlieren method. Fluid motions can also be “visible” to radio 
waves: an article by Ottersten (1969) displays remarkable photographs of atmospheric 
vortices imaged by radar in clear air. These include bouyant vortex rings, roll cells aligned 
with the wind in a stably stratified layer, and Kelvin-Helmholtz rollers in the tropopause. 
Perhaps someday we will set? radar photographs of aircraft wakes as well. 
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An important practical problem is that the vortex wakes of large aircraft pose a haz- 
ard to following aircraft and therefore, during instrument take-offs and landings, fixed 
distances are maintained between aircraft. These spacings are thought to be too conser- 
vative most of the time and to increase airport capacity, NASA lias set itself the goal of 
developing a system which dynamically sets the spacing (Perry et al. 1997). Several ele- 
ments comprise this system including computer simulation of wakes and vortex detection 
from the ground. 

In the NASA effort, two methods of vortex detection are being concentrated upon: 
doppler lidar at infrared frquencies and radar. Lidar is the more mature of the two 
efforts and relies on scattering by aerosols such as v/ater droplets which may begin to 
form in the low pressure of the vortices even when cloud is not visiblef. Field tests at 
airports using 10.6 /mn lidar (Campbell et al. 1997, Kopp 1994) have detected vortices 
at ranges of up to 300 m; in clear weather this limit is set by the depth of focus of the 
optics. Hannon & Thomson (1994) report ranges of up to 4 km at 2.09 fim using high 
energy pulsed lidar. 

There are two motivations for considering radar. First is the concern that lidar cannot 
work in fog and rain due to increased absorption. Radar, on the other hand, is unaffected 
by rain and fog at sufficiently low frequencies. The second factor is that while lidar has 
adequate range to protect the modestly sized approach corridor currently set by the 
planners, radar has a potentially greater range. 

1.2. Previous theoretical work 

The results of radar tests that have detected aircraft wakes are all interpreted using the 
theory (Tatarski 1961) of scattering by refractive index fluctuations caused by atmo- 
spheric motions in clear air. The theory is based on the fact that a turbulent velocity 
field in the presence of mean vertical gradients of potential temperature and humidity 
lead to fluctuations in refractive index. The velocity field is assumed to be statistically 
homogeneous and isotropic locally and to lie in the Kolmogorov inertial range. Refractive 
index gradients are therefore also homogeneous and isotropic. One result of the analysis 
is that if the wave number of the incident radio wave is fc, the only wavevector of the 
refractive index capable of scattering to the observer is: 

k = k(\ — o), (1.1) 

where i and o are the unit vectors in the incident wave direction and the direction from 
the target to the observer Respectively. This is the so-called Bragg condition. In particular 
for backscattering (6 = — i) one gets: 

k = 2A;?. (1.2) 

Thus, the only flow wavevector which backscatters is aligned with the radar beam and 
has twice the radar wavenumber. According to Tatarski’s theory the radar cross-section 
per unit volume of isotropic turbulence in the inertial range is (e.g. see Ottersten 1969) 

t) = 0.38O- 1/3 , (1.3) 

f While scattering is the usually proffered explanation, it is interesting to note that according 
to Carlon (1970) water aerosols are capable of strong absorption and emission in the 8-13 /i 
range and that “water droplets in this wavelength region behave much more as absorbers and 
emitters and much less as Mie scatterers that has been imagined.’’ Two of three vortex lidar 
tests were 10.6/a and therefore re-emission could very well be present. This fact suggests that 
passive radiometry of condensed water aerosol in the vortex cores might be another method of 
detection. 
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where A is the radar wavelength and r; has units of L ~ { . The radar cross-section can 
thus be obtained if one knows C„, the so-called second-order structure constant of the 
index of refraction fluctuations. Equation (1.3) has become widely quoted in the aircraft 
vortex detection context as well, with the belief that aircraft wake turbulence creates a 
turbulent energy cascade similar to atmospheric turbulence. Equation (1.3) is also used 
to reduce experimental data (Chadwick et al., Nespor et at. 1994). In particular, the raw 
radar cross-section (RCS) reported by the radar is divided by the volume of the pulse to 
obtain q and then the value of calculated using (1.3), is plotted. 

Marshall & Myers (1996) attempted to predict C “ as a function of position in a 
Reynolds averaged calculation of a statistically 2D vortex wake by using an algebraic 
closure relation for C 3 in the inertial range in terms of gradients of mean velocity and 
potential refractive index. This model is analogous to that given by Ottersten (1969, p. 
1184) for atmospheric turbulence. 

There are also classified reports (with unclassified titles) cited by Gilson (1992) to 
which we were not privy. 


1.3. Previous experiments 

(i) Noonkester k Richter (1980) detected the wakes of aircraft taking off (f = 3 GHz, 
A = 10 cm, range R =100-300 m). Their FM-CW radar, in contrast to a pulse radar, 
transmits and receives continuously allowing for a small minimum range. An interesting 
sinusoidal pattern was observed in the trace of target height (above a fixed ground 
location) vs. time. They interpreted this to be hot engine exhaust spiraling around each 
vortex. Whether engine exhaust might be put to use for vortex detection during approach 
is briefly discussed in §6.2. 

(ii) Using a similar radar, Chadwick et a/. (1984) detected the wakes of both departing 
and arriving aircraft but do not distinguish between the two in their results. The radar 
used was rather modest (200 W of power, 8 ft diameter dish) and the ranges were R ^ 1 
km. The distribution functions of RCS in the atmosphere ahead of the plane and the 
wake behind the plane overlapped a little and for positive identification of the wake, the 
authors suggested use of the spread in doppler velocity. They suggested two mechanisms 
for scattering. First, they argued that aircraft vortices created a turbulent cascade similar 
to those assumed in the theory of atmospheric scattering (Tatarski 1961) but with a C\ 
higher than that of the ambient atmosphere. They found no directional sensitivity which 
lent some credibility to the the assumption of isotropy. Secondly, they suggested that 
heat and moisture from the engine further increases reflectivity. 

(iii) Using a 1 MW pulse doppler radar (/ = 5.6 GHz, R = 2.7 km), Nespor et al. (1994) 
detected the vortices of a small fighter aircraft in approach configuration looking axially 
along the vortex. 

(iv) So far, Gilson (1992) has the best documented and controlled experiment. Us- 
ing powerful pulse doppler radars having 2-7 MW of peak power, he measured the radar 
cross-sections of the wake of a C-5 A aircraft at a range of 15 km. Returns were detectable 
at five frequencies between 0.162 GHz and 5.7 GHz but not at 35 GHz. From the observa- 
tion that the radar cross-section was relatively flat as a function of frequency and dropped 
off somewhere above / = 5.7 GHz, he concluded that particulates were not involved (they 
would give f 4 Rayleigh scattering). He noted that the frequency dependence was not the 
Kolmogorov Z 1 / 3 . The RCS decreased with altitude 8a id vanished above the tropopause 
from which he concluded that the scattering mechanism was related to “low altitude 
climate.” Fluid mechanicians will anticipate that the relevant fact here is the decrease 
with height of the Briint-Vaisala frequency, the appropriate measure of stratification. He 
found that the RCS (measured 1 km behind the plane) was insensitive to engine thrust 
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ami the RCS (measured at 10,000 altitude and 0.6 km behind the plane) was insensitive 
to flap setting. He thus concluded that the main mechanism was “turbulent mixing by 
the wake vortices of existing atmospheric index of refraction gradients.*’ 

(v) For completeness, we mention a test (lanrmzzelli et al. 1998) in which a radar was 
used to detect index of refraction variations caused, remarkably, by sound waves launched 
as a tracer into the aircraft wake from the ground. The sound waves are scattered by the 
vortex flow and the doppler shift they induce on a radar beam can be measured. This 
technique is an extension of a routinely used method for wind profiling and is known as 
RASS (Radio Acoustic Sounding System). 


1.4. Present work 

Led by the observation that turbulence is suppressed in an aircraft vortex, §2 considers 
three simple laminar and largely two dimensional mechanisms for producing refractive 
index variations. 

Given the index of refraction, §3 performs a scattering analysis using the Born weak 
scattering approximation. It is argued that the usual assumptions of a uniform incident 
wave and a far-field are not generally valid for a long and coherent target such as an 
aircraft wake. Therefore the shape of the incident beam and full Green’s function are 
taken into account. This analysis has to be implemented numerically and so to provide a 
simple analytical result, and for insight into the results, a simplified analysis is presented 
in §4. This analysis makes the approximations of a far-field and a cylindrical beam that 
is wider than the cross-sectional width of the target. 

Section 5 presents numerical predictions of the radar cross-section for the 2D mech- 
anisms, comparing the approximate analysis, the full Born analysis and experiment. 
Section 6 presents some practical reccomendations including a sample design for an ST 
type radar to detect the vortices themselves. 


2. Three mechanisms for inhomogeneities in dielectric constant 

The three sub-subsections which follow discuss the change in dielectric constant due to 
three effects: (i) radial density gradient in the core of each vortex, (ii) three dimensional 
instabilities in the vortex cores, and, (iii) transport by the vortices of the atmospheric 
density. Following Gilson (1994) we assume that particulates are unimportant to the 
scattering mechanism in clear air. 

The refractive index n of humid air for frequencies below 20 GHz is given bv Thaver 
(1974): 

(n - 1) x 10 6 = 77.6 + 64.8 + 3.776 x 10 5 . (2.1) 

Here p a is the partial pressure (mb) of dry air, p v is the partial pressure (mb) of water 
vapor and T is the temperature (K). The first two terms are due to the induced polar- 
ization of air and water molecules, respectively. The third term is due to the permanent 
dipole moment of the water vapor molecule. p v can be obtained from the relative hu- 
midity, RH = Pv/Psat, where p sat is the saturation pressure of water vapor and can be 
obtained as a function of temperature using the formula in Reid et al. (1987, p. 757). The 
dielectric constant is simply s = n 2 . 
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2.1. Radial density yradient in the vortex cores 

The radial momentum equation (for laminar flow) together with the assumption of ho- 
mentropic flow can he integrated to yield the radial density variation in each vortex: 


P{r) 

Poo 


r™ 1 V \ t^t 
l-(7-l)/ -±dr) . 

Jr r c oo / 


( 2 . 2 ) 


Here and p 0 G are the ambient speed of sound and density, respectively. Equation (2.2) 
applies to both the density of water vapor and the density of air. 

Two profiles for the circumferential velocity, V# were considered. The first is that of 
the Lamb-Oseen vortex: 




(2.3) 


The following values, provided by Rennich (1997) as representative of a large commercial 
aircraft were used: core radius, r 0 = 4.724 m and circulation, T = 526 m 2 s _1 . We take 
the core radius, ro, to be the location where the tangential velocity has a peak and for 
the Lamb-Oseen profile we have 6 = r 0 /l. 12141. With these values Eqs. (2.2) and (2.3) 
give a small density dip: p(0)//?oo = .998. 

The second profile of circumferential velocity was provided to us by Dr. P. Spalart of 
Boeing and represents his fit to flight observations: 


Here 


Vo 


r 

27T7}bo 


1189 t) 2 , for Tj < .0103, 

£(1.27 4- log(^)/4) -14 + lj , otherwise. 


(2.4) 


,£ v < 2 - 5 > 

We chose bo = 47.88 m, the value of vortex spacing provided by Rennich (1997). This 
profile gives a greater minimum: p(0)/poo = 0.988. Figure 1 plots the relative change 

= — - 1 (2.6) 


in dielectric constant from the ambient value (e 0 ) for the two velocity profiles and for 
two cases of ambient relative humidity (dry and 80%). The ambient temperature and 
pressure were taken to be 15° C and 1013.25 mb, respectively. We-observe that humidity 
has a weak influence and the profile of Spalart has a more pronounced dip of Ae r at the 
vortex axis « -6 x 10 -6 . 


2.2. Three dimensional fluctuations in the vortex cores 

Little is known about turbulence within an aircraft trailing vortex but it has come to 
be appreciated that both the cascade to small scales and turbulence production are 
drastically reduced by rotation. Turbulence models currently used in aerodynamic as 
well as atmospheric codes are unable to capture these effects. 

With respect to the first effect, the rate of energy transfer and hence dissipation is 
reduced (see, for example, Bardina et a/. 1985) because rotation scrambles non-linear 
interactions. 

A more important effect is the drastic alteration of turbulence production. There are 
three mean velocity gradients which can act as turbulence production mechanisms within 
the vortex: 

(i) There is a jet/wake flow' along the vortex axis but this decays quite rapidly behind 
the wing. 
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Figure 1. Radial distribution of Ae r inside a vortex. : Lamb-Oseen vortex, (RH)^ = 0%; 

: Lamb-Oseen vortex, (RH)oo = 80%; : Spalart’s profile, (Eq. 2.4), (RH)oc = 0%; 

: Spal art’s profile, (Eq. 2.4), (RH)oo = 80%. 

(ii) There is a circumferential velocity around the axis. In the central part of the core 
the production of turbulent kinetic energy is zero because of solid body rotation. As one 
moves away from the center, fluid elements are sheared as they rotate. The Bradshaw 
criterion for stability of shear in a rotating frame when applied to a vortex with single 
sign vorticity predicts that the flow is locally stable at each r (LeBlanc & Cambon 1997). 
Indeed Qin (1998) found via direct numerical simulation of turbulence within the vortex 
that there is negative production of turbulent kinetic energy. 

(iii) Finally, the only remaining mechanism for turbulence production is the strain 
induced by one vortex upon the other. It is responsible for the commonly observed Crow 
(1970) instability which has a wavelength of about nine times the vortex separation. A 
shorter wavelength instability, called the Widnall instability in the context of vortex rings 
is also possible at sufficiently high Reynolds numbers. This instability has been observed 
in the laboratory by Thomas & Auerbach (1994) and Leweke & Williamson (1998). 

The question arises whether even shorter waves can be excited by the strain in the 
presence of the stabilizing effect of the rotating shear. To address this question without 
performing a full stability analysis of a vortex pair we consider a stability analysis which 
accounts for only the local velocity gradient. Such an analysis is valid for waves much 
smaller than the vortex core size. The local velocity gradient following a fluid element 
has three parts: rotation, strain induced by the other vortex and rotating shear. The first 
two alone lead to the elliptic instability (e.g. Landman & Saffman 1987). To this we add 
rotating shear. The velocity gradient matrix then becomes: 

( 0 —7 — e 0 \ / — sin fit cos fit cos 2 fit 0 \ 

7 — e 0 0 1 4 - 5 [ - sin 2 sin fit cos fit 0 ]. (2 7 ) 

0 0 0 / V 0 0 0 ) 

The first term represents the elliptic streamline flow which has rotation and strain while 
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Initial azimuth of wavevector ((f) 

Figure 2. Amplification factor of disturbances at various locations in a Lamb-Oseen vortex 

subject to strain by another vortex : r/6 = 0.0; : r/S = 0.01; : r/S = 0.10; 

: r/6 = 0.20; : r/6 = 0 30. 


the second term represents a shear (of rate 5) rotating with angular velocity ft. The 
value of fi was chosen to be the rate of rotation of fluid particles in the elliptic flow, 
i.e. ft = (7 2 - e 2 ) 1 / 2 . Had ft been chosen to be the rate of rotation, 7, of fluid particles 
without the strain then the shear axes would get ahead of the axes of a fluid element 
which is not physical. For an axisymmetric circumferential flow «#(r) we have 


7 = — and 5 
r 


( du 9 

V ^ 



The strain rate due to the other member of a vortex pair is 


( 2 . 8 ) 


e — 2 x 
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(2.9) 


where the induced strain has been multiplied by 2 in order that the growth -ate and 
ellipticity of the elliptic streamline flow match that of a vortex of uniform vorticity sur- 
rounded by potential flow. The evolution of a single mode was considered at several values 
of r/6 inside a Lamb-Oseen vortex with a value of 6 /bo = .088 given by Rennich (1997) 
for a large aircraft in approach. The modes investigated were those known to be the most 
unstable for the elliptic streamline flow. These modes have wavevector lying in a cone 
at 60° relative to the vortex axis. Within this cone the initial azimuthal angle y? is free. 
Figure 2 shows the amplification factor of the velocity perturbation during the time that 
the vortex pair descends five vortex spacings. The different curves are for different loca- 
tions, r/S within the core. At the center of the core (chain dashed line) there is no shear 
and the elliptic instability is active for all c p. As one moves outward, the shear stabilizes 
all modes except those located in a thin band near p — 90°. This wavevector direction 
corresponds wavefronts that are parallel to the shear flow and therefore not deformed 
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by it. By r/6 = 0.3 (long-dashed curve) the growth of even these waves is suppressed 
by the shear. We thus conclude that roughly only the inner 20% of the vortex core, or a 
diameter of 170 cm for a large airplane, is capable of supporting an elliptic instability, 
and even this produces small amplification factors. The instability will be subject to a 
viscous cut-off. Using the analysis ofLandman & Saffrrian (1987) one estimates that for a 
large airplane the cut-off occurs for scales < 14 cm. Down in scale from the cut-off there 
may be smaller scales due to the non-linear cascade. 


2.3. Dielectric constant variation due to atmospheric stratification 

In the literature it has been postulated that the scattering mechanism for aircraft vortices 
is turbulent mixing of atmospheric index of refraction gradients 7 ’ (Gilson 1994), the 
turbulence being assumed to be similar to atmospheric turbulence. This sub-section tries 
to develop a concrete picture for the quoted phrase. 

2.3.1. The 2D picture without velocity fluctuations 

Suppose that the vortex pair descends through a distance Az = z f -z from its altitude 
of formation, z/. As it descends it carries, within an oval, ambient atmospheric fluid from 
z /. The oval has semi-axes 1.0456 o and 0.8656 0 (Milne-Thomson 1968, p. 360) which we 
shall approximate in the scattering analysis by a cylinder of circular cross-section of 
radius a = 0.95b o . 

Let tildes denote fluid properties within the oval and let overbars denote properties in 
the ambient atmosphere. Then introducing the total pressure, p = p a + p V) into (2.1), 
converting to SI units and eliminating pressure in favor of density using the ideal gas law 
we get the following equation for the difference in index of refraction between the oval 
and the surrounding atmosphere: 


[n(z) 


Pv(z) 


■n(z)]xl0 6 = 223 [p(z) -T>(z)\- 36.7 [£„(*) -p v (z)\+ 1.08 x 10 6 

T(z) T(z) _ 

As it descends, the fluid in the oval compresses adiabatically in response to increasing 
ambient pressure. A density difference, p(z) - p(z ), from the ambient value will remain 
only if the ambient atmospheric density gradient is different from the adiabatic one. This 
difference is usually expressed in terms of the parameter N, the Brunt Vaisala frequency 
defined so that for small Az (Pedlosky 1979) 


= (211) 

The first term in (2.10), which represents the dry air contribution, is thus determined. 

Since the pressure of the descending fluid is equalized to the ambient we can obtain 
its temperature using Eq. (2.11): 


t-t(z) = +TeTLl AZ ' 


(2.12) 


Tbobtam p v (z), the ambient vapor density, we introduce the ambient relative humidity 
RH(z) as an input parameter so that p v (z) can be obtained from 


am _ Pyjf) _ rh(z) Pmi (t) 

P(z) P{z) P( z ) ’ (2 ' 13 ) 

where the function Psnt (T) is given in Reid et a/.(1987, p. 657). Finally, to obtain p v (z), 
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Case 

RH/ 

RH(z) 

N ins' 1 

•Ae, 

1 

0 

0 

0 

0 

2 

0 

0 

.017 

— !).8 x IO" 7 

3 

.60 

.60 

0 

-2.8 x 10~ 7 

4 

.60 

.60 

.017 

-i.i x io -6 

5 

.592 

.603 

0 

-4.4 x io -6 

6 

.592 

.603 

.017 

-2.8 x 10 -6 


Table 1. Dielectric constant in developed in the oval for various atmospheric conditions after 

1.5 r units. 


use the fact that all constituents of the descending air compress by the same factor so 
that: 


Py ( Z f ) Pv( z ) 

P( z f) P( z ) ’ 

Then introducing the relative humidity at formation 
obtains: 


(2.14) 


(RH /) as an input parameter one 


Pv(z) 


RH fp 9J a(T f )p(z) 
P(z) 



( p(z) 

\P(z) 9 9 



(2.15) 


The value of N near sea level varies from average summer levels of 0.014 s' 1 to mild 
winter values of 0.02 s 1 to a value of 0.03 s' 1 for a Fairbanks winter (W.D. Kriese, 
Boeing Aircraft, Private Communication; Pedlosky 1979, p. 331). Thus, winter weather 
produces greater reflectivity than tropical weather (for the same vortex descent altitude 
and neglecting humidity). What is the range of Az that is relevant for the wake hazard 
problem? The lifetime of trailing vortices has been measured to range between 1.5 and 
8 r units (Spalart 1998, his fig. 1) where one r unit corresponds to the time takes for a 
vortex pair to descend one vortex separation distance. Hence the minimum separation in 
time between aircraft that could ever be contemplated (in the absence of techniques to 
destroy vortices) is 1.5 r units. A vortex separation distance of bo = 47.88 m implies that 
the smallest detectable voitex descent should be 72 m. This assumes that stratification 
neither impedes nor speeds the descent. 

Table 2.3.1 shows Ae r % 2An for various _atmospheric conditions taking Az = 61 m, a 
slightly more conservative value. We used F(z), p(z), and T(z) values for the Standard 
Atmosphere at sea level. Case 1 is the trivial case (neutral and dry armosphere). Case 2 
is for stably stratified dry air. Case 3 has a neutral atmosphere with uniform humidity. 
Case 4 combines the stratification of case 2 with the humidity of case 3. We see that the 
total effect is not quite additive. Case 5 uses the vapor pressure for the model atmosphere 
given in the Handbook of Geophysics (1960, p. 13-5) in which the air is slightly drier at 
the higher altitude. The resulting A£> is higher than for uniform humidity. Case 6 has 
both stratification and humidity gradient. 


2.3.2. Dielectric constant of the oval in Gilson s experiment 

" To secure a value for the dielectric constant Ae r for the experiment, performed on 
Kwajalein atoll, we used the density and index of refraction vs. height sounding given by 
Gilson. Since gradients relative to adiabatic are needed, one quantity had to be assumed 
to obtain them. This was the temperature at flight altitude (5000 ft) which we took to be 
295 K assuming a balmy sea level temperature of 85° F and a lapse rate of - .005 K/m. 
By comparing the measured index of refraction with the value for dry air at the local 
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density we gathered that 7% relative humidity must have been present at flight altitude. 
The value of the Brunt-Vaisala frequency N was determined to be .017 s _l by comparing 
the actual density gradient with the adiabatic one. We need the descent altitude of the 
ova! at the first measurement instant 66 secs after aircraft passage. Using the flight speed 
and nominal aircraft weight provided together with the theory of an elliptically loaded 
wing, the vortex circulation and vortex separation were inferred to be T = 387.6 nr/s 
and bo — 53.4 m. For vortex pair descent unaffected by stratification, this yields a descent 
of 77 m in 66 secs. The parameter determining the effect of stratification on vortex pair 
descent is N 0 = 2iri\bo/r (Spalart 1996). Its value is 0.79 and using the simulation 
results of Spalart (1996, his figs. 4 and 10a) we concluded that that stratification would 
shorten descent from 77 m to about 64.6 rn. Thus we obtained Ac r = -8.4 x 1CT 7 . 

2.3.3. The boundary of the oval 

(i) Laminar boundary Thermal diffusion smooths the gradient in Ap at the boundary 
of the oval and the flow sweeps away and replenishes the diffused Ap. This may be locally 
modelled as a strained diffusion layer: at equilibrium we have (e.g. Leonard et al 1987): 

Ap(r?) = Ap(l) erf(7?/<$ diff ), S dlff = (2/c/e) 1/2 (2.16) 

where rj is normal to the interface, k is the thermal diffusivity and e is the strain rate. 
For a Prandtl number of 0.7 and 

6 = 2tt(6o/2) 2 ’ (2 ' 17) 

which represents the strain rate at a distance of half the vortex spacing, one gets 

/T \ ~ x/ - 

<5di ff = 2.1&0 ( - 1 , (2.18) 

which equals 2.2 cm for Gilson’s experiment. 

When flaps are deflected (as during a landing) multiple vortices are created on each 
side of the wake and their interactions can lead to periodic ejection of fluid from the oval 
and engulfment of ambient fluid. This leads to a stretching and folding of diffusion layer 
as shown in the simulations of Spalart (1996). However the minimum thickness of the 
folded layers cannot be smaller than the estimate (2.18). 

Any vorticity fluctuations within the core will induce potential flow fluctuations in the 
outer flow. This will also cause the boundary of the oval to become distorted. Maxworthy 
(1970) provides experimental evidence for this in the context of vortex rings and Rom- 
Kedar et a/.(1990) provide a geometric description. 

(ii) Turbulent boundary In the presence of stratification baroclinic torque generates 
a shear layer at the boundary of the oval (Spalart 1996). Can this shear layer become 
turbulent? To answer this, the Reynolds number Re* = ( AU/2)S/i / was estimated and 
compared with a critical Reynolds number. Here A U and 6 are the velocity jump and 
laminar thickness respectively of the shear layer and were estimated as follows: (i) Spalart 
(1996) provides a model (his Eq. 6), calibrated using simulations, for the circulation 
generated in the shear layer. This allows one to infer the average circulation per unit 
length (hence the average velocity jump A U) along the oval perimeter, (ii) The shear 
layer thickness, set by bala nce between viscous diffusion and vortex induced strain, is 
given by S ss n l / 2 b 0 / y/TJu. 

The Reynolds number was found to grow linearly with descent height and even with a 
descent a quarter of what it is at 66 secs, we have Re s = 73 for Gilson’s (1992) experiment. 
According to the analysis of Betchov k Szewczyk (1963) for a shear layer, this is quite 
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unstable with an e-folding time of about 2 secs. At this rate, transition to turbulence in 
the shear layer would take about 20 secs, hi reality, the shear layer is curved. However, 
since the baroclinic circulation is opposite that of the primary vortices, the circulation 
decreases outward and so curvature can only further destabilize the layer. Thus we are 
convinced that the boundary of the oval is turbulent in Gilson’s measurements. 

Hie circulation in the Kelvin-Helrnholtz vortices can be estimated as the circulation 
in one wavelength of the instability. From it the strain rate in the braid region can be 
obtained and thence the thickness of the braids. At 66 secs, <5 b raids = 0.45 cm. This is 
smaller than S^ff we expect it to determine the the wavelength at which the RCS cuts 
off. 

2.3.4. Evidence for the role of stratification in radar reflectivity 

Here we discuss facts immediately accessible from the experiments. Later in the paper 
we will perform a more quantitative comparison. 

(i) Gilson observed increasing RCS with distance behind the plane. Nespor et al.( 1994, 
p. 658) noted that the cross-section tended to increase as the vortex dropped into the 
lower beams. Their respective interpretations were developing turbulence and growth 
of the turbulent region (so that a greater fraction of the radar pulse volume consisted 
of turbulence). However, the observations are more simply explained by the increasing 
refractive index contrast between the oval and the ambient. Indeed RCS oc (Ae r ) 2 cx 
(Az) 2 and (Az) x t (when the descent is not impeded by stratification). Thus, the radar 
cross-section should increase by a factor of 4 (6 db) for every doubling of the distance 
behind the plane. Several plots in Gilson (1992) show an increase that are a little smaller 
than this value probably because the descent is impeded by stratification. 

(ii) Gilson notes that the RCS decreases with altitude. This is consistent with the 
decrease in the Briint-Vaisala frequency, A f , with altitude. 

In conclusion, there is experimental support that atmospheric stratification is playing 
a role in experiments in the simple way we have described. Since stratification has sea- 
sonal, diurnal and geographic variations its use as a basis for detection would have to be 
combined with atmospheric measurements and would be unusable in neutral conditions 
or when the vortex descent is strongly impeded by stratification. 


3. Scattering Analysis with the Born Approximation 

3.1. General formulation 

The scattering analysis for an aircraft wake has two peculiar features that; ought to be 
recognized. (1) The wake is a long target which always fills the incident beam in one 
direction and is quite coherent along this direction. Thus one cannot follow the type of 
analysis performed for a collection of rain drops, for instance, where it is assumed that 
at each point within the beam, a plane wave of uniform intensity and polarization is 
incident on each drop and the scattered powers of all the drops simply add because the 
scattered fields are uncorrelated. In the present situation, one has to account for the 
variation in beam properties across the whole target. (2) More troublesome is the fact 
that the usual far-field assumption cannot be made because as the distance from the 
radar to the target is made larger, the illuminated size of the target also increases! More 
concretely, the condition that should be satisfied in the far or Fraunhoffer zone is: 




A R r 


«: 


l, 


(3.1) 
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where R r is the distance from the receiver to the target, L s is a characteristic size of 
the scattering region and A is the wavelength. The difficulty in satisfying (3.1) for a 
beam filling target arises from the fact that L„ = Rt A# / 2 grows with range Rt from 
the transmitter to the target. Here A# is the half-power beamwidth. In particular, for 
co-located transmitter and receiver (3.1) becomes: 


7r/? r (A/9)- 

etr = — A — - « 1. 


(3.2) 


Ironically, the far-field condition becomes more difficult to satisfy the further one is 
from the target. The experiments of Chadwick et a/.(1984) and Gilson (1992) we have 
6 ^ Sj 15. Tatarski’s (1961) analysis makes the assumptions of a far-field and a 
uniform incident field. Hence it is not formally valid for a coherent wake. Fortunately, 
the problem remains tractable even if these assumptions are relaxed. 

VVe consider time harmonic electric and magnetic fields 


E(x, t ) = Re [E(x)e -iu,t ] , H(x, t ) = Re [H(x) e - iL “'] , (3.3) 

where E and H are complex. Rationalized MRS units ( Jackson 1962, p. 611) will be 
used to write the electrodynamic equations. The starting point is the solution to the 
scattering problem as an integral equation (Ishimaru 1978, vol. 1, p. 16): 


Here 


n 5 (x)= / Ae r (x')E(x')G(x,x')dx\ 

Jv 

. . .. e(x') 

Ae r (x ) = — — - - 1. 


£o 


(3.4) 

(3.5) 


e ifcA 

G(X,X) = 4^A’ A = I x " x 'I> (3-6) 

is the free space Green’s function and II S is the Hertz vector of the scattered field: 


E s — V x V x n s , H s = -iw£ 0 V x n s . (3.7) 

To calculate E and H from II requires evaluating curls (with respect to x) of the 
integrand in (3.4). We list them here for completeness. Denoting by C = Ae r (x')E(x') 
the vector in the integrand which is constant with respect to x we get: 

V x CG = VG x C, ( 3. 8) 

V x V x CG = V x (VG x C) = (C • V)VG - CV 2 G. (3.9) 

Equations (3.8) and (3.9) involve first and second derivatives of G which are: 
dG 1 e ifcA 

dii “ 4jt A 2 Xi ~ 1//A ) ’ ( 3 - 10 ) 

d 2 G 1 e‘* A r/ ., 

diidxj = + lkx i x J A ~ 2x i x j) (>* - 1/A) + xnj/A] , (3.11) 

where 5jj is the Kronecker delta. 


3.2. The Born approximation 

The integrand in (3.4) contains the total electric field (incident + scattered) which is in 
general unknown, however, since in the present case Ae r = O (l(T 6 ) one can invoke the 
Born approximation that it can be set equal to the incident field. This corresponds to 
retaining only the first term in a Neumann series solution to the integral equation. The 
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y (thumb) 



Figure 3. Wake coordinate system. Use the indicated fingers of the left hand to follow the 
transformations (3. 15)— (3.19) discussed in the paper. 


Born approximation is valid when the scattered electric field is much smaller than the 
incident field in the region of the target. Equations (3.4), (3.9) and (3.11) imply that the 
following three conditions have to be satisfied for this to be true: 

Ae r { 1 , kL s , (kL s ) 2 } < 1 . (3.12) 

It should be noted that Ishimaru (1978) gives only the second condition perhaps because 
he was thinking of scattering by particles smaller than a wavelength. In that case ArL s < 1 
and the second condition is the more stringent one. 

3.3. Calculation of the radar cross-section 

The definition of radar cross-section RCS (e.g. Skolnik 1970, p. 1-4) can be written as: 

RCS — Air ft 2 P° wer flux density spatially averaged over the receiving aperture 

r spatial maximum of power flux density across incident beam at target ‘ 

(3.13) 

Power flux density has dimensions of power per unit area. The numerator in (3.13) will 
be evaluated at a single point which is valid if the size of the receiving aperture is smaller 
than R ry an assumption which is difficult to violate. Both numerator and denominator 
will be evaluated as time averages. The time averaged scattered power flux density vector 
is 

Ss(x) = ^ Re (E s (x) x H s *(x)} , (3.14) 

where Re denotes the real part. The numerator in (3.13) was evaluated as |S 9 (x)|. As a 
check, the numerator was also evaluated as Si(x) • o where o is a unit vector (defined 
in (3.25) below) from t he target to the receiver. A difference between the two was never 
detected. 

Appendix A. 2 shows that the denominator in (3.13) is simply (l/2)ar 0 - 
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3.4. Coordinate systems 
We choose the coordinate system of x and x' in the integral (3.4) to be located in the 
wake with a convenient origin. In particular for an axisyrnmetric target the origin is 
placed at the target axis as shown in figure 3, its axial location being chosen to coincide 
with the axial location where the beam centerline intersects the xy plane. In particular 
the beam centerline is chosen to intersect the xy plane at (x = 0,7 / = t/o), //o being the 
beam pointing error. 

Appendix A obtains the electric field in the incident beam given the electric field on 
a surface (the aperture) near the antenna. The incident field is given using coordinates 
centered on the aperture. In the body of the paper we use (x*, y m , z m ) for aperture coor- 
dinates. The plane of the aperture is x*y * while 2* points normal to the aperture. We 
now relate beam coordinates (x*, y* , z*) to wake coordinates (x y y,z) by performing a 
series of rotations and translations. These transformations serve to define the viewing 
angles (a*, ft) of the target from the transmitter. Similarly, (a r ,ft) are viewing angles 
from the receiver. 

Imagine that the transmitter is looking at the wake from the same side as the reader is 
looking at figure 3. This figure also shows which fingers of the left hand the reader may 
employ to visualize the transformations. The reader’s palm will eventually become the 
plane of the aperture, the middle finger will become the direction of propagation, and 
the first finger will become the direction of polarization. The first transformation turns z 
(the middle finger) away from the reader so that it can eventually become the beam axis. 
The x-axis (first finger) is also simultaneously flipped. The second transformation shifts 
the origin to (x = 0, y = yo) to ensure that the beam axis intersects the xy plane where 
we want it. The coordinates resulting from the first two transformations are therefore 



Next, rotate the coordinate system by angle ft € [0,2 tt] around x^ 1 ) axis. This is a 
counter clockwise rotation about the first finger. Thus ft is an elevation angle to the 
target. We have 


/ x< 2 > 

v* (2) 


i o o 

0 cos ft — sin ft 
0 sin ft cos ft 



(3.16) 


Next is a rotation by angle a t € [ — tt/2, zr/2] about the axis. This is a counter- 
clockwise rotation about the thumb and we get 


/ x< 3 > 

y l3) 

\ z< 3 > 


cos at 

0 

-sin a t \ 

( x< 2 > 

0 

1 

° 

w (2) 

sin at 

0 

cos a t ) 

V 


(3.17) 


Next, translate the coordinate system in the negative z^ direction by the range R t from 
the transmitter to the point x = (0,7/o?0) on the target: 


x ^ 
j/ (4) 

zW 




(3.18) 


Finally, to the get the polarization vector oriented correctly, the coordinate system is 
rotated by angle 7 about the z^ axis (about the middle finger, counter-clockwise looking 
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into the palm): 


( cos 7 sin 7 0 

— sin 7 cos 7 0 

0 0 1 


(3-19) 


The x M axis (first finger) now points along the polarization direction. Putting all the 
transformations together gives: 

x* — cos 7 {— xcosa* — sin a t [(y - yo) sin fit — z cos fi t ]} + 

sin 7 {(y - y 0 )cosft + z sin fit) , (3.20) 

y m = cos 7 {(y - y 0 ) cos fit 4- z sin fit} - 

sin 7 { — xcos a* — sin a* [— z cos fit + (y — yo) sin fi t ]} , (3.21) 

z* = Rt — x sin ct t 4- cos at [—z cos fit 4- (y — yo) sin #*] . (3.22) 

The incident electric field is in the x* direction and to obtain its (x, y, z) components for 
use in the integral (3.4) we note that 


x* = — cos a t cos 7X 4* (cos fit sin 7 - cos 7 sin a* sin fi t ) y 4* 
(sin fit sin 7 4- cos 7 sin at cos fit ) z, 


(3.23) 


from (3.20). Just as the “incident direction” i from the transmitter to the point (0,yo,0) 
is 

i = z* = — sin QjX 4- cos a t sin fity — cos a t cos fi t z, (3.24) 

we define the “observer direction” o from the point (0,yo,0) to the receiver in a similar 
way but using angles a r and fi r and reversing the sign of the vector: 


o = sin a r x — cos a r sin fi r y -I- cos a r cos fi r z, 


(3.25) 


The location of the receiver is then: 


x r = yo y 4- RrO. 


3.5. Numerical integration procedure 

The numerical integrations required to evaluate the integral (3.4) were performed using 
two methods. 

Method A: This is the more general but less efficient of the two methods. It uses 
Simpson’s rule in cylindrical coordinates (x, r, <fi). The integral over the infinite direction 
x was truncated at |x| = 3«b where Ub = Rt sin#b Is the characteristic beam radius at the 
target, #b being the half beam width. In the radial direction, the integrals were* truncated 
at r = 3 S e when the vortex was the target. When the oval is the target, the integrand 
is zero for r > a. The number of quadrature points in each direction was determined 
from the input parameter N q i which is the number of quadrature intervals desired per 
characteristic length / of fluctuations of the integrand. The characteristic length in each 
direction was chosen as follows: 


l x = min (A, u b ), 1$ = min(A, tt b ), l r = rnin(A, u b , S r ) y (3.27) 


where 


f (S f , if the vortex is the target, 
\ (l/2)a, if the oval is the target. 


(3.28) 


If the length of an integration is L then Nt = int (L/l) 4- 1 characteristic lengths cover 
the integration length and the number of quadrature intervals was determined as = 
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i\ f} t max(3, /V/), the purpose of the max being to ensure a minimum number of quadrature 
points regardless. For each calculation of radar cross-section presented for tfiis method 
we began with ~ 6 and kept on doubling it until two successive doublings showed no 
perceptible change in the radar cross-section. Computational cost prevented use of this 
method above / ^ 0.16 GHz. 

Method B: This method was developed to efficiently calculate the radar cross-section 
for the oval at the relatively high frequencies of Gilson’s (1992) experiment. The imple- 
mentation is currently specialized to backscattering for the oval (and for the oval with a 
slice of ambient fluid; see figure 8). For brevity we will describe the method for normal 
incidence and 3 t = 0. 

Each of the six integrals needed to evaluate the scattered E and H fields is of the form: 


I = J Vi z)e l0(x ' y ' z) dxdydz , (3.29) 

where z is the direction along the beam axis and the cross-beam coordinates are (x,y). 
The basic idea of the method is to divide the integration interval in a given direction 
(say C) into many sub-intervals and to approximate the amplitude and phase in each 
sub-interval as a linear function of 


where 


•'Ci Jo 


(3.30) 


A.4 = .4(6) - .4(0), = <j>( C->) - 0(Ci). 


The integral in each sub-interval can then be obtained analytically: 


(3.31) 


A <fi 


AC 


,4(0) + -A.4 


^(Ci 


4* 


A Kl ) + -A A 


p i«HCa) 




(3.32) 


The integral along the beam ( z direction) was done first and tests revealed that one 
sub-interval was sufficient to obtain an accurate integral for each region. For the sliced 
oval, one interval for each piece was used. The resulting integral over an interval has 
the form (3.32) with £ replaced by z. Now the two terms in the result (3.32), when 
considered as a function of the remaining coordinates, are again in the basic form (3.30) 
with amplitudes and phases a function of x and y. Hence, the same method can be 
applied successively in these directions. In these directions hundreds of sub-intervals 
were required for convergence. 


4. Approximate scattering analysis 

4.1. General formulation 

To obtain analytical insight into the numerical results and as a check on the numerics for 
special cases, it is useful to have closed form expressions which approximate the radar 
cross-section. Such expressions are developed in this section and, later in the paper, their 
predictions will be compared with the formulation of the previous section. 

Four assumptions will be invoked. The first assumption is that of the far-field which 
holds when (3.1) is satisfied. In this case the scattered electric field implied by (3.4) is 
given by (see Ishimaru 1978, vol. 1, p. 17): 

E s (x) = f(o,?)e i ** p /fl r , 

- Ar f 

f(o,i) = — Jv {-8 x [8 x E(x / )]} &£ r (x')e- ikx '- B dx\ 


(4.1) 

(4.2) 
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The second assumption is the Born approximation of the previous section: E(x') — > 
Ej(x') in the integrand. 

The third approximation replaces the conical incident beam by a cylindrical beam in 
the vicinity of the target as explained in Appendix A. 2. The incident electric field then 
has the form 

Ej(x) = x'F(u;«r)e'***, (4.3) 

where 

u = (x' 1 + y’ 2 ) / , (4.4) 

is the normal distance from the beam axis. Substituting for from (3.20) and using 
(3.24), the incident electric field takes the form 


Ej(x) =x m F(u]t 7)e lkxi e lk *, 


where 


(4.5) 


4> = R t -y 0 cos a* sin /?* . (4.6) 

Two types of beam shapes, justified in the appendix, will be considered. The first type 
is Gaussian: 

F( U ;<T)=e-“ a /- 3 j (4.7) 

and corresponds to a Gaussian illumination of the transmitting aperture. The second 
type corresponds to that produced by a uniformly illuminated circular aperture: 


F(u\<t) = 


2 Mu/ a) 

(u/<t) 


(4.8) 


Using (3.20) we have 

u 2 =i‘ 2 + y* 2 = Ax 2 + Bx(y - yo) + Cxz + D(y - y 0 ) 2 + E(y - y 0 )z + Fz 2 , (4.9) 


in which 

A = cos 2 a f , B — sin 2a* sin/?*, C = — cos /?* sin 2a* , 

D = cos 2 fit + sin 2 a* sin 2 /?*, E = cos 2 a* sin 2/?*, F = cos 2 /?* sin 2 a* + sin 2 /?*. 

(4T0) 

Keeping ail the terms in (4.9) makes the resulting integrals analytically intractable. Thus 
we make the fourth approximation that the beam is much wider than the cross-sectional 
dimensions of the wake. In the wide beam limit, the beam amplitude varies only along 
the axis of the wake and is constant within each cross-section of the wake. Thus setting 
y = z = 0 in (4.9) we get. 


u 2 = x* 2 + y ml = Ax 2 - Bxy 0 + Dyl . 


(4.11) 


Substituting (4.5) into (4.2) gives: 


where 


f(S,i) = 


sA:“e' 




47T 


f ±e r (xL')F{u;a)e lkx ' {Us) 

Jv 


(4.12) 


s = -o x [3 x x*]. (4.13) 

We have |s| = sin \ where \ is the angle between x* and 6. 

To calculate the scattered power flux we need the scattered magnetic field correspond- 
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ing to (4.1). Using the Maxwell equation 


one obtains 


H s ' = — — V x E,, 
lUJflQ 


(4.14) 


H * = 


_ — i kR r 


Rr^Ho 


(3 x f)’ 


(4.15) 


Substituting (4.15) into (3.14) and noting that according to (4.12) f is of the form f = sC 
where s is a real vector perpendicular to o and C is a complex number, we conclude that 


S s (x) 


2/?;ca/io 


|f| 2 S. 


(4.16) 


The incident flux in the definition of the radar cross-section (3.13) is taken to be that 
for the original conical beam, namely, ce o/2. Hence (3.13) gives simply: 


RCS = 47r|f|-\ (4.17) 

Finally, let us keep ready at hand the following geometric quantity which appears in 
Eq. (4.12): 

x • (i - o) = -Px - Qy - Sz, (4.18) 

where 


Begin by Fourier 
tions: 


P = sin a t 4- sin a r , 

(4.19) 

Q = ~ cos at sin fit — cos a r sin fi r , 

(4.20) 

S = cos a t cos fit + cos a r cos fi r . 

(4.21) 

4.2. Scattering off a Fourier mode 



decomposing the dielectric constant in the axial and azimuthal direc- 


Af r (x) = Ae r (r;A: I ,m)e i( ** I+m * ) . (4.22) 

To justify the sum in k x instead of an integral, pretend that there is a periodicity length 
L x and then let L x oo. 

Substituting the Fourier expansion (4.22) and the Gaussian incident beam (4.7) into 
(4.12) gives: 




fc x ^ T71 


where 


sk 2 e* k<tk 

f(k x ,m) = — — — I x rdr Ae r (r; k x ,m)I^{r), 

and / 0 and f x are the following integrals (for a Gaussian beam): 

f 1 * 

[$= I dip e~ lkr ( s cos 

Jo 

fO © 

I x =e~ D Vo/ cr / dx Q~( Ax2 - Bx v*)I a2 Q-' xx i kp - k *) 
J — oo 


(4.23) 


(4.24) 


(4.25) 

(4.26) 


To perform the integral with respect to <f> let kQ = fcsin^t, kS = kcos<t>i, so that 
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k = k(S- + Q 2 ) 1 / 2 and 

r‘2n _ 

U = J d<t> e~'* r ) . (4 27) 

Using an integral tabulated in Gradshteyn & Ryzhik (1980) (p. 482) we get: 

h = 2e‘ m0, i m 7r./ m (-*r). (4.28) 

To evaluate /*, let A: = A:/ 1 — A: x and complete the square on the exponent to get 


j jo 

I r = ?—jL e (yo/<rf[B : ‘/(4A)-D]-ik , Byo/C2A)-{k'vf/(4A) 
A 1 / 2 ^ 


(4.29) 


Note that there is a Gaussian cut-off factor if the beam is off-center: exp{[fl 2 /(4.4) - 
D \(Vo/v) 2 h where the argument can be shown to be negative definite. This corresponds 
to the ‘blip” as the radar scans across the vortex. Setting y 0 — 0 one obtains 


^ c L‘2 ikR t 1/2 _ 

f(k x ,m) = — 

' 2 A 1 /' 1 


• m f 

Jo 


rdr Ae r (r; k x ,m)J m (—kr) (4.30) 


This has a cut-off due to finite beam size as represented by the factor e ~( k ' <r /( 4/ U 
(A: = kP—k z ). It implies that the contribution to the radar cross-section versus frequency 
by a fluctuation in refractive index at wavenumber k x is a “bump” which is centered at 
radio wavenumber k = k x /P and which has a half- width of Sk = 2A 1 / 2 /(Pa). This result 
is analogous to the Bragg condition of Tatarski (1961). The cut-off is eliminated when 
it has infinite width, i.e. when P = sina r + sin a* = 0. This holds if the receiver is at 
the direction of specular reflection (a r = -a t ). Note that for backscattering (a r = a f ), 
perfectly normal incidence (a r = a t = 0) is required to eliminate the cut-off. The cut-off 
is also eliminated if a r = a t ± jr and it corresponds to the receiver being on the other 
side of the transmitter relative to the wake. 

For the Bessel function beam the analysis is very similar; only the integral I x is differ- 
ent: 


I, = fl d* = *1 

J - OO A 


[\-{k'a) 2 /A] ir2 , 0 ^ k'a/A 1 / 2 < 1, 
0, k'a/A ‘/ 2 ^ 1, 


(4.31) 


where u — (Ax 2 ) 1 / 2 . In deriving (4.31) the beam offset y 0 was set. to zero for analytical 
tractability. Thus one sees that the Gaussian cut-off for the Gaussian beam has been 

replaced by a sharp cut-off and the same conditions hold for the elimination of this 
cut-off. 


4.3. Approximate scattering by three mechanisms 

Each of the following three sub-sub-sections specialize the approximate scattering anal- 
ysis to each of the three mechanisms discussed in §2. 

4.3.1. Radial density gradient tn the vortex 

We will represent the radial variation ( k x = m = 0) of dielectric constant shown in 
figure 1 by a Gaussian: 

Ac r (r) = Agfexp (-r 2 /S‘“). 

Then the radial integral in (4.30) is 


(4.32) 
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This has a cut-off due to finite vortex core size which is eliminated when k = k(S~ -f 
Q ) ^ — 0. This is satisfied for three cases. The first is for incident and scattering direc- 
tions that are along the vortex. We eliminate this as a possibility since the assumption 
that the beam strength does not vary in the cross-section of the target becomes invalid. 
The second case corresponds to a t = a r and ,3 r = 8 t + tt. The third case corresponds to 
forward scattering (a r = ~a t , fl r = + tt). For this case, the cut-off due to beam width 

is also eliminated and therefore this is the optimal arrangement. 

4.3.2. Three-dimensional fluctuations within the vortex 

As discussed in §2.2 three-dimensional fluctuations should arise in the inner part of the 
vortex cores from the elliptic instability. The azimuthal mode excited by this instability 
is rn = 1 and for this mode we model index of refraction using known forms for velocity 
eigenfunctions: 

Ae T (r;k x ,m = 1) = eA^J^AvrJe - ' 3 ^ 3 , (4.34) 

where e is the amplitude of the disturbance and k r is the radial wavenumber which we 
set to k r = y/3 k x . This value corresponds to the unstable wavevectors which lie on a 60° 
cone. The resulting radial integral in Eq. (4.30) is given in Gradshteyn & Ryzhik (1980) 
(p. 718, 6.633.2): 


[ r = -€A£?6~ex p 



(4.35) 


For large arguments the modified Bessel function I\ has an exponential factor and there- 
fore the last two factors in (4.35) can be combined to give a Gaussian bump with respect 
to frequency: 




exp 


yj nk r k 

The center and half width of the bump are given by 


f (* r -*) 2 


(4.36) 


V3k x 

(S 2 + Q 2 ) l /2’ 


A* = 


2 /S e 

(5 2 + g2)i/a- 


(4.37) 


We also have a Gaussian bump in Eq. (4.30) due to finite beam width. In order that the 
two bumps not annihilate each other they must be very close. In particular the reflectivity 
is maximized when they occur at the same frequency. The condition for this is: 


S 2 + Q~ = 3P 2 (4.38) 

For backscattering Eq. (4.38) simplifies to cota* — 3 or a = ±30°. Thus the radar 

beam must be aligned with the wavevector of the instability. As P 0, i.e. as we 
approach specular reflection the peak of the beamwidth bump moves to higher frequency 
and its width becomes larger. In order to maintain reflectivity the second bump (Eq. 4.36) 
must also move out to higher frequency. This is accomplished by letting (S 2 + Q 2 ) 1 / 2 0, 

i.e. by approaching forward scattering, The increase in frequency yields an increase in 
radar cross-section due to the leading k 2 factor in (4.30). 


4.3.3. Oval 

As mentioned in §2.3.1 the oval can be approximated a s a cylinder of radius a (A: x = 
m = 0). In the absence of fluctuations and thermal diffusion the refractive index will be 
constant inside the cylinder and zero outside. In this case the radial integral in Eq. (4.30) 
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IS 


f (t ~ n \p° 

Ae" / rdr./ 0 (Ar) = 

■A) 


(4.39) 

1 he integral is intractable if we wish to treat error function diffusion layers. Following 
(4.33), however, we estimate the cut-off due to both laminar and turbulent layers at the 
boundary of the oval by the following factor applied to (4.39): 


-*V/ 4 


(4.40) 


where 6 is either S,utr or Irakis (§2.3.3). 


5- Numerical Results 

The next three sub-sections contain numerical results for the radar cross-sections of 
the 2D mechanisms discussed in §2. Since at this time amplitudes and wavenumbers of 
3D disturbances in the vortex cores are not known, plots for this mechanism are omitted. 

5.1. Scattering by 2D density variations in a vortex 

The results of the full Bom analysis were obtained using Method A (§3.5) and a Gaussian 
beam having a half- width of 0^ = 1-5° and angle of polarization 7 = 0 is assumed. 

In (4.32) the width S e of the well in dielectric constant was taken to be .0256o and the 
depth As® of the well was taken to be -6 x 10“ 6 (both from figure 1). 

Figure 4a plots the radar cross-section at a range of 1 km. The results of the approx- 
imate analysis (lines) agree quite well with those of the Born analysis (symbols) even 
though the far-field parameter becomes 0(1) (Figure 46). The general features of the 
curves are in accord with the discussion in §4.3.1. One observes (i) a low frequency cut-off 
for back-scattering away from normal incidence due to finite beam width (chain-dashed), 
(ii) a higher frequency cut-off at the scale of the vortex for back-scattering at normal 
incidence (dashed) and specular reflection (dotted), and (iii) elimination of both cut-offs 
for forward scattering (solid). The long-dashed curve shows that some increase in fre- 
quency and reflectivity can be obtained even somewhat away from the forward direction, 
in particular when the transmitter and receiver are both looking up at the vortex with 
an elevation of 30°. 

The limit frequency for the Born approximation is obtained by setting L s = max(u b , 6 e ) 
in (3.12). The result is / <$: .74 GHz. This is only marginally satisfied for the rightmost 
portion of the figure. From the Mie result for scattering by a sphere we expect that if 
the Born approximation were not made, the main difference would be that the solid line, 
which increases linearly, would begin to oscillate and decrease in average' slope as the 
weak scattering approximation began to be violated. 

To emphasize differences between the full Born and approximate anatysis the range 
was varied beyond 1 km. For a point target the radar cross-section should be independent 
of range but not so for an extended target. For very short ranges the vortex will fill the 
beam and we expect the RCS to increase quadratically with range. When the range is 
large enough that the vortex fills the beam only in one direction then we expect that the 
radar-cross should increase linearly. To verify this, the frequency was kept fixed at the 
peak for normal back-scattering (/ = 5 x 10 7 Hz) and the range was varied (see circles in 
figure 5a). For off-normal back-scattering the frequency was fixed at / = 2 x 10 6 Hz. The 
square symbols in figure 5a show that as the beam radius increases with range the cut-off 
becomes sharper. It is interesting to note that the dependence on range according to the 
approximate analysis (dashed line) is quadratic, that is, as if the target were beam filling 
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Figure 4. (a) Radar cross-section for a vortex due to a Gaussian beam at a range of 1 kin. 
The symbols correspond to the Born analysis while the lines give the result of the approximate 
analysis. , □ : back-scattering away from normal incidence (a t = a r = 30°, 0 t = 0 r =0); 

, o : back-scattering at normal incidence (a t = a r = = j5 r = 0) o * 

S PTl"o re £ ection (“< = -«r = 30°, fit = 0r = 0); , v : nearly forward scattering 

with 60 offset in 0 (a, = -a, = 30°, fi t = O,0 r = 120°). , a : forward scattering 

(art — a r - 30 ,/? { - 0,^ r = 180°). (b) Far-field parameter, ex for the top figure. 


in both directions. Figure 56 shows that the approximate begins to show deviations when 
the far-field parameter becomes greater than unity. 

5.2. Scattering by the oval 

5.2.1. Comparison with Gilson's (1992) Experiments 

Experimental procedures and parameters (Gilson 1992, Roth et al 1989) are now dis- 
cussed: 
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FIGURE 5. (a) Dependence of the radar cross-section of a vortex on range. The symbol- 
s correspond to the analysis with only the Born approximation while the broken lines give 

the result of the approximate analysis. , o : back-scattering at normal incidence 

(a* = a r = 0 yflt = 0 r = 0), / = 5 x 10 7 Hz; , a : back-scattering away from nor- 
mal incidence (a t = a r = 30°, ft = A = 0), / = 2 x 10 6 ; , reference lines to indicate 

quadratic and linear dependence, (b) Far-field parameter, ex for the top figure. 


(i) Gilson (1992) reports RCS data in db m (i) 2 per meter of wake length (Gilson 1992, 
his p. 37) along the radar beam. Hence, we multiplied his RCS values by our estimate 
(101.5 m) for the diameter of the oval. This value is consistent with the size of the wake 
in the radar returns. 

(ii) The range of the target was 15 km and we chose the case having the lowest flight 
altitude (5000 ft) presented in Gilson (1992). The radar was pointed normal to the wake. 

(iii) All the radars in the experiment transmit a right circularly polarized (RCP) wave 
and receive both right and left circular polarized (LCP) waves. So far our analysis has 
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Band 

A0 


fo 

/it in 

Waveform 

Pulse width, 

FM spread, 


(rnrads) 

(degrees) 

(GHz) 

(GHz) 

ID 

r (/xs) 

A/ fm (MHz) 

VHF 

48.8 

2.37 

0.162 

0.084 

V.25C 

0.25 

0 

UHF 

19.2 

0.934 

0.422 

0.21 

u.tc 

0.10 

0 

L 

10.6 

0.516 

1.32 

0.39 

L2 

2 

20 

S 

5.2 

0.253 

2.95 

0.79 

S3 

3 

60 

c 

5.2 

0.253 

5.67 

0.79 

NB 

10.2 

6 

K u 

0.76 

0.037 

35 

1.0 

N2 

50 

12 


Table 2. Gilson’s radar parameters. AO: half-power beamwidth; d h : e" 1 half-width of a 
Gaussian beam; /i im : Limit frequency for validity of the Born approximation. 



Figure 6. Comparison with Gilson’s (1992) experiment, o , Gilson’s data at 66 secs behind the 

plane, • : Full Born analysis for oval; * , : Full Born analysis for oval with a slice 

of entrained fluid. ■ : Approximate analysis for the oval. 


been for linear polarization. However, a CP wave is simply a superposition of two linearly 
polarized waves 


E(x) = £(x)(e! ±ie 2 ), (5.1) 

H(x) = H (x) (ei ± ie 2 ) , (5.2) 

where e* J_ e 2 . First note that the power flux of a CP wave is simply twice the power 
flux of one of its linearly polarized components. If the incident beam has the form (5.2) 
then the linearly polarized analysis can be done separately for each component and the 
resulting scattered fields added. In this manner, it was verified that for all the cases 
presented here, the resulting scattered wave is also CP (but of opposite helicity). Thus 
both incident and scattered powers are doubled compared to the linearly polarized case 
and the RCS is left unchanged. 

(iv) Gilson used a pulsed radar whereas our analysis is for a time-harmonic wave. 
Appendix B describes the treatment of a pulse and Table 5.2.1 gives the characteristics 
of the waveforms. 
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Relative range, c{t — Rr/c) /2, in meters 

Figure 7. Returned pulse envelope calculated for Gilson’s (1992) VHF case. The abcissa is 
time expressed as a range, c (t - R r /c) /2; the division by 2 converts from round-trip travel 

distance between scattering centers to one-way distance. : Oval; : Oval with slice 

of entrained ambient fluid. 


Figure 6 compares predicted and measured RCS values. The open circles are measure- 
ments 66 secs behind the plane. The predictions of the full Born analysis for the oval are 
the solid circles. Each value plotted represents the peak value of RCS(f) in the returned 
pulse. The predictions agree with the measurements for bands L, S, and C. 

At VHF and UHF, the pulse widths are so short that in the calculation, the two 
scattering centers of the target (see discussion after Eq. B 14 in Appendix B) produce 
two distinct echoes (see solid line in figure 7). If the electric field of the two echoes 
are added, or if a long pulse is used such that the two echoes overlap, then the result 
matches the experiment. Such addition can only be justified if signal processing in the 
experiment added returns from multiple range gates within the wake but according to W. 
Gilson (Private communication) this was not so. Addition could also occur if there are, 
in reality, more than two scattering centers in the target spaced sufficiently close to each 
other. To test this suggestion, a slice of width S ent = .20 a was removed from the oval to 
represent entrained ambient fluid (figure 8). With elevation angle 0 = 5.8°, the value for 
Gilson’s (1992) experiment, the resulting echo, shown as the dotted line im figure 7 now 
consists of one main pulse. The RCS predictions are shown as the solid triangles in figure 
6. The agreement at VHF has now improved but not at UHF which has a very narrow 
pulse width. It is possible that the target still has more internal structure we have not 
accounted for. 

Gilson obtained no return at K a band (35 GHz). Figure 9 shows the attenuation in 
RCS, due to the oval boundary estimated using (4.40). For the thickness of the turbulent 
layer the result (solid line) fits observations: there is very little attenuation at 5.66 GHz 
but a lot at 35 GHz. The thicker laminar layer, however, attenuates too strongly (dotted 
line). 

An upper frequency limit fy Hn for validity of the Born approximation was defined using 
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Figure 8. Model for the oval with a slice of entrained ambient fluid. 



Figure 9. RCS attenuation due to the oval boundary. : 5 = (J diff = 2.2 cm; ■ 

= ^braids = 0.45 cm; o : Frequencies corresponding to Gilson’s experiment. 


(3.12) by replacing the inequality with an equality and using 

L s = ma x(R t sin(0 b ), a) (5.3) 

for the characteristic radius of the scattering region. Table 5.2.1 shows that Gilson’s 
experiment has frequencies between factors of 2 and 6 above /, im . Therefore in the future 
it would be worthwhile to consider another approach such as geometric optics which valid 
for large ka but arbitrary Ae r . 

The analysis predicts a rapid cut-off away from normal incidence. On the other hand, 
comparing experiments at 45 and 90° incidence (Gilson’s runs 1 and 17, respectively) 
does not reveal significant differences. Moreover, Gilson searched for a flash at normal 




Radar reflectivity of aircraft vortex wakes 27 



Figure 10 . Test of various approximations. : full Bom analysis; : with the far-field 

Green’s function; : with the wide cylindrical beam approximation; : with both ap- 
proximations. - : same as uie wide cylindrical beam approximation but without assumption 

(d) in Appendix A. 2. The thicK lines depict linear and quadratic dependence on range. 


incidence and found none. We estimated in §2.3.3 that the boundary of the oval is very 
likely to be turbulent and this may account for the lack of directional dependence in the 
experiment. 

5.2.2. Test of various approximations 

Range dependence. The effect of various approximations was studied using Method 
B by providing the option to use either the far-field or the lull Green’s function and the 
option to use the either the conical beam or the approximate cylindrical beam. Figure 
10 shows the results. The parameters chosen were those of Gilson’s S-band test except 
that the range was varied. The result (solid line) of the full Bom analysis shows the R 2 
behavior for small ranges, characteristic of a beam filling target, followed bv the linear 
behavior characteristic of a target that fills the beam in only one direction. As the range 
increases both the far-field G (dotted) and approximate beam (dashed) results seem to 
agree with each other and to overpredict the RCS of the full Bom analysis (solid). This 
agreement partially understandable since both require smallness of a far-field parameter 
(see approximation (d) in Appendix A. 2) in order to be valid. If we drop assumption 
(d), the remaining assumptions ((a)-(c) in Appendix A.2) should begin to hold as the 
range increases and the long dashed curve bears this out. When both the far-field G and 
approximate beam are used, the result (chain-dashed) gives an R 2 behavior throughout 
which is erroneous for large R. This curve agrees, as it should, with the result of the 
approximate analysis; the comparison is not shown for clarity. 

Frequency dependence. Figure 1 1 shows the frequency dependence of the RCS of the 
oval for a time harmonic incident wave. The rest of the parameters are those for Gilson’s 
V HF case and his data (circles) has been included for reference even though they are 
strictly for a pulse. The solid line is the result of the full-Born analysis and shows a flat 
frequency dependence. The curve is punctuated by oscillations. These oscillations have 
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Figure 11. Frequency dependence of the RCS for the oval. : full Born analysis (normal 

incidence); : approximate analysis (normal incidence); : far-field G (normal inci- 
dence); : full Born analysis (15° off-normal incidence); : approximate analysis (15° 

off-normal incidence); o : Gilson’s experiment. 


a regular period (versus frequency) of Af = c/(4a) due to interference between the two 
effective scattering centers of the target (see Appendix B). Except for the lowest frequency 
decades, the number of plot points is so coarse that the fluctuations appear irregular. 
Using the far-field G produces the dotted curve which only slightly overpredicts the RCS. 
The curve for the wide cylindrical beam is very close to the dotted curve and is not shown 
for clarity. However, when both far-field G and wide cylindrical beam are invoked (to 
yield the approximate analysis) the result (long-dashed curve) has an erroneous RCS oc k 
behavior. For the same number of plot points, the dotted and long-dashed curves have 
the same pattern of oscillations indicating that the period and phase of the oscillations 
is the same; only the amplitude is higher. The other two curves with the sharp cut-off at 
lower-left are for off-normal incidence (a r = a* = 15°). 


6. Practical Recommendations 

6.1. Prospects for detection of the vortex cores 
No test so far has been at frequencies as low as 49 MHz where, according to the anal- 
ysis, the density gradient in the vortex cores reflects (see dashed line in figure 4). Since 
this mechanism is independent of atmospheric conditions, it is worthwhile to discuss its 
potential use in air traffic control (ATC). 

(i) It is fortunate that 49 MHz is where many radars operate around the world for 
atmospheric profiling. Such radars are called ST or MST radars and they employ an 
antenna array. Thus the predicted reflectivity could very well be tested. 

(ii) Rain clutter is not a problem at this frequency, though clear air turbulence might 
be. 

(iii) Most ST radars are used for measurements in the high atmosphere but recently 
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various problems associated with measurements at lower heights have been overcome and 
a boundary layer profiler has been developed (Vincent et ai 1998) which has a height 
coverage of 400 - 3200 rn. This radar achieves a range resolution of 150 m with a pulse 
width of l fin. 

(iv) For ATC use, the required cross beam resolution (bearnwidth) should be a fraction 
of the minimum airplane spacing expected (about 2 nautical miles). The smaller the 
bearnwidth the larger the required array but, since the antenna gain is larger, less power 
is required. Desiring both modest array size and modest power we impose a cross beam 
resolution of 0.3 n.m. at ranges < 3 km resulting in a 10 x 10 array with elements spaced 
3 m apart and a power of only 360 W. Since this value is surprisingly modest, details of 
the calculation follow. 

For an antenna array the half-power bearnwidth (in degrees) broadside to the array is 
(Skolnik 1970, p. 11-2): 





( 6 . 1 ) 


where N e \ is the number of elements, and the gain on transmit and receive is: G t = G r « 
irN e \. The transmitter power required for detection is: 


W^VpxSNR 

1 G t G r A 2 L sys x RCS’ 1 ' 

where jV p = —185 db VV and L sys = -7 db are the noise power and system loss for an 
ST radar (Kingsley & Quegan 1992, p. 199). SNR = 20 db is the desired signal to noise 
ratio, A = 6 m is the wavelength and RCS = —77 db, a value we calculated for a vortex 
at 3 km using the bearnwidth of the 10 x 10 array. 

(v) ST radars have limited steering capability, usually just a few fixed beam directions 
but this may be sufficient to observe the approach path of a few adjacent runways at 
a fixed distance from the touch-down point. Since it is unlikely that three dimensional 
disturbances will exist in the vortex across a broad spectrum of wave angles, detection 
may be possible at only normal incidence. 

(vi) Finally let us note that the oval will also scatter at 49 MHz. 

The above estimates have ignored smaller aircraft types and detection of vortices of air- 
caft in landing connfigurations. Such considerations would have to be taken into account 
in a full-fledged feasibility study. 


6.2. A note on the role of engine exhaust 

At approach, the exhaust of a large aircraft has a temperature of 477 K and each engine 
emits 1200-6000 Ib/hr of water vapor (P. Spalart, Private communication). If we consider 
the temperature contribution alone (the first term in Eq. 2.1) we get Ae r = -2.1 x 10“ 4 . 
This is roughly two orders of magnitude greater than for the density variation in each 
vortex and for atmospheric density stratification and would lead to a 40 db greater radar 
cross-section. However, Gilson (1992) found that the wake RCS at 1 km behind the plane 
did not change when the engines were run at idle or full power. This is probably due 
to turbulent mixing of the exhaust with the ambient. Diversion of some engine exhaust 
into the laminarising flow of the vortex cores may allow the exhaust to maintain its heat 
further downstream. 

A simpler alternative is to detect the exhaust itself as in the Noonkester & Richter 
(1980) experiment. Using the radar trace of the exhaust spiral (their figure 5) we deduced 
that the circulation of the vortices is 400 nr/s. In clear weather, an even simpler alter- 
native is to view the exhaust spiral in infrared in order to deduce the vortex circulation. 
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To estimate the contribution of water vapor in the exhaust, use the upper value of 
6000 Ib/hr/engine and assume an exhaust diameter of 2 m. This yields a partial pressure 
of vapor equal to I . I mb and a contribution of \e T of 3.7 x 10 -6 which is much smaller 
than that due to temperature. 


7 . Conclusions 

(i) The measured radar cross-sections in the experiments of Gilson generally agree with 
the values caclulated using the simple picture of a oval carrying atmospheric fluid as it 
descends. The lack of directional dependence is not predicted by the analysis, however, 
estimates indicate that this is because the boundary of the oval is very likely turbulent. 

(ii) The density gradient in the vortices themselves produces index of refraction vari- 
ations comparable to those produced by the oval. This provides a scattering mechanism 
independent of atmopsheric consitions and could be detected using an ST type radar at 
about 49 MHz. 

(iii) The use of a far-fieid Green’s function is not a priori valid for a wake that is 
coherent along its length and therefore a more general analysis was developed. The use 
of a far-field G was found to produce a constant db error in the radar cross-section 
for large ranges. The neglect of beam conicality produces similar errors. When both 
assumptions are made, the dependence of radar cross-section on range is altered from 
the correct linear behavior (as expected for a target that fills the beam in one direction) 
to quadratic behavior. Similarly, the frequency dependence is altered from being flat to 
cx k. 

We thank Drs. W. Gilson and P. Ingwersen (MIT Lincoln Labs.) for providing us with 
information on the experiment, Dr. P. Spalart (Boeing) for suggesting that we account 
for density variation inside a vortex core and for supplying the information referred to 
in the text. Thanks are also due to Prof. J. Jimenez (Univ. Madrid) for suggesting the 
inclusion of water vapor in the refractive index and to Dr. M. Wang (CTR) for help with 
integrals. 


Appendix A. The electric field in the incident beam: Ej (x, y, z) 

The incident electric field is obtained using aperture antenna analysis (Silver 1949) 
where one derives the far diffraction field from knowing the field on a surface (the “aper- 
ture”) near the antenna. We shall take the aperture to be circular. 

For convenience we shall, in the manner of a subroutine, use a dummy coordinate 
system (x, y, z ) centered on the aperture. The plane of the aperture is xy and 2 is normal 
to the aperture. These coordinates should not be confused with those in the body of 
the paper where they refer to wake coordinates and where (x* , y* , z m ) refer to aperture 
coordinates. 


A.l. Conical beam 

According to Silver (1949), the electric field created in the far zone by a plane aperture 
A on which the electric field is specified to be E = ££(x')x is the following integral over 
the aperture (Silver 1949, p. 167, his Eq. 124): 

Ei(x) = -^-xe ifcr (l + cos6>) f £“(x')e i * x '- x/l * l cix'. 

4?r r J A 


(Al) 
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Here (r, <9) are the spherical polar coordinates of x, 9 being the polar angle measured 
from the ^ axis. Equation (A l) is derived starting from vector Kirchhoff integrals and 
introducing several approximations valid in the limit of small wavelength compared to 
aperture diameter in which case the significant portion of the energy is confined to small 
9. In addition, the phase of the field on the aperture has been assumed to be constant. 

For an axially symmetric distribution of E'f on the aperture the integral in (A 1) 
becomes 


r A = 2irj gdg E z (g)Jo(kg’ sin#), (A 2) 

where g‘ = yj x' 1 -f i/-. Two aperture illuminations are considered. For the case of a 
uniform illumination: 


we get 


( Eq j Q ^ Qa > 

[ 0, otherwise, 


(A3) 


I A = 2*Eoe n — 7 — 7 —. 

kg a sin v 

For the case of Gaussian illumination of the aperture 

E n z {g) = E 0 e-^^\ 


(A 4) 
(A 5) 


we get 

I A = vcr'iEo exp (-^/c 2 c7 2 sin2 » (A 6) 

thanks to an integral tabulated in Gradshteyn & Ryzhik (1980, p. 717, 6.63.4) To calculate 
the radar cross-section according to (3.13) requires the maximum value at the target of 
the time averaged power flux Si = l/2Re (Ej x H'}, where H, can be obtained from the 
Maxwell equation H; = 1 /(iw^o)V x Ei. If in evaluating VxEj from (A 1) we retain 
only the leading order term ( O (1/r)) and introduce = l/(c 2 e 0 ) we get 


Si = 


o c£ o {E' x ) 


2 Z ^ 

“Z, 

r 


(A 7) 


where E l x denotes the x component of Ej. The radar cross-section is independent of 
the amplitude Eq since both incident and scattered powers are proportional to Efi. For 
convenience, we chose 


( 2R t /(ihg*), for uniform aperture illumination, 

\ 2 R t /(ik(rl), for Gaussian aperture illumination, 


(A 8) 


in order to make equal to e ikRt the maximum value of E\ with respect to 9 at fixed 
r = R t . The denominator in (3.13) is then 


maximum power density in incident beam at target = -c£ 0 - 


(A 9) 


In the calculations, instead of specifying the aperture radius as the input parameter, we 
specify the half bearnwidth 9 b and from it calculate the implied value of g a . For the 
Gaussian beam 9^ is defined to be where I A drops to e“ l of its peak value and for the 
Bessel beam it is defined to be the location of the first zero thus giving 


3.83171 _ 2 

Arsing’ ° a Arsing 


(A 10) 
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A. 2. Approximate cylindrical beam 

The cylindrical beam for the approximate scattering analysis in §4 of the paper is obtained 
after introducing four approximations in (A 1): (a) Replace cos# by 1. This is valid for 
a small angular beamwidth. (b) Replace r in the denominator by a constant range Rt . 
(c) Similarly, in (A 2) replace sin# = u/r by u/R t , where u = \J x 1 + y 1 is the normal 
distance from the beam axis, (d) Replace e lkr with e lk= . Approximations (b)-(c) are valid 
when the relative change in beam width across the near and far ranges of the target is 
small. For approximation (d), however, a condition similar to the the far-field condition 
(3.1) must hold: 


irL; 

xt^ L 

Introducing the stated approximations into (A l) gives: 


(All) 


E = 2^r xe '''" (A12 > 

Substituting (A 4) and (A 6) for I ^ into (A 12), replacing sin# with u/R t , using (A 10), 
defining u^/Rt = sin#b, and substituting values of Eq from (A 8) we get 


E = 


xe 


ikz 


2 Ji(u/a) 
u/cr 


a = ^b/3.8371, 


for the Bessel function beam and 


(A 13) 
(A 14) 


E = xe tkz exp (-u 2 /a 2 ), (A 15) 

a = u b , (A 16) 

for the Gaussian beam. For the incident power flux in the cylindrical beam we use the 
same value, given by (A 9) as in the original conical beam. 


Appendix B. Scattering by a pulse 

Since the scattering problem is linear, treatment of a pulse is simply a matter of Fourier 
superposition of the results of the time-harmonic analysis. The analysis was performed 
for back-scattering at normal incidence. 

We begin by reducing the originally linear vector problem to a linear scalar problem. 
This will provide the set-up for using formulae already derived in Ishimaru (1978). 

For the time harmonic analysis the incident electric field (Eq. A 1) has the form: 

Ej(x, f; w) =x Re [Fi(x, u;)e~ lu/ *] . (B 1) 

As in Appendix A we are using x to denote aperture coordinates. For the accompanying 
magnetic field we assume the same spatial dependence: 

Hi(x, t;u) = c£ 0 y Re [Fi(x, uj)e~ iuJt ] . (B2) 

According to Maxwell’s equations E is related to the curl of H and therefore E and 
H should in general have different spatial dependence. However, to the extent that the 
wavelength is much smaller than the beamwidth, derivatives of F(x, cj) normal to the 
beam axis will be much smaller than derivatives along the beam axis and the forms (B 1) 
and (B2) will nearly satisfy Maxwell’s equations and be divergence free. 

Henceforth we will write the incident field as a function of time only, taking it to be 
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evaluated at the point x = (0 ,t/o,0) where the beam centerline intersects the centerplane 
of the target. By a Fourier synthesis of (B 1) and (B‘2) in time an arbitrary L-> time 
dependence can be constructed: 


Ei {t) = x Re [ui(f)e- iu, °'] , (B3) 

Hj(£) = cs 0 y Re [u i (*)e“ lu ’ 0 '] , (B4) 

wliere we have chosen to write the time dependence as the product of a carrier e~ tuJot 
and a complex modulation u*(£). The complex modulation is itself written as a product: 

«»(*) = /(t)e* {t) , (B 5) 

where f(t) and 4>(t) are real and describe amplitude and phase modulation, respectively. 

We now assume that for the time harmonic analysis the scattered field evaluated at 
the receiver is also of the form (B 1) and (B2) with replaced by F s . This assumption 
was verified to be true from the numerical results. Therefore by Fourier superposition we 
can write for arbitrary time dependence: 


E s(t) = xRe [u s (t)e“ iu;o *] , (B6) 

H s 0) = ce 0 yRe [u s (t)e" iu, ° f ] , (B7) 

Our incident and scattered fields thus depend only on the scalar functions ui(t) and ti s (£) 
and we are thus set-up to use Ishimaru’s (1978, vol. 1, p. 94) discussion for scattering by 
a pulse. Using the fact that scattering is a linear and causal problem, Ishimaru derives 
the following “input/output” relation for the complex modulation of the scattered pulse 
in terms of the modulation of the incident pulse: 



H (u) 4- u)q , t) e loJ * did. 


(B 8) 


where 

1 f°° 

U 'M = ^ J Ui(£)e ,u7 * dt. (B 9) 

has the interpretation that if u\{t) = 1 (i.e. we have the time harmonic case) 
then u s (t) — H(uj 0} t). In the time harmonic analysis we chose the incident field to have 
a normalization such that Ui(t) = e lkRt . Hence H((jjQ,t) is simply the u 3 we calculate in 
the time harmonic analysis divided by e**^*. Since we have considered a steady target, 
thus losing the doppler effect, H(oj 0 ,t) does not depend on time. 

Let us obtain H(ui) implied by the result of the approximate analysis for the oval. For 
backscattering at normal incidence and using the asymptotic form of the Bessel function 
for large ha , Eqs. (4.30) and (4.39) give: 


+ /?,.) / o \ 

E s (x) = x — — (ka) [ / 2 <j&£® cos f 2*a - -n J (B 10) 

The bandwidth of the pulse is typically small compared to the center frequency and so 
ka can be considered to be constant when outside the cosine. Hence, over the bandwidth 
of the pulse 


where 


H(uj) = Ae ikR " cos 



k = u>/c, 


(B 11) 


A£?<7(Ar 0 a) 1/2 

W r 


(B 12) 
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is a real constant independent of u and k$ is the wavenumber corresponding to the 
center frequency. For the hill Born analysis of the oval, it was found that H (uj) had a 
behavior approximated by (B 11) but with a different (and complex) value of .4 which 
was obtained by fitting the form (B 11) to the numerically obtained results. 

Substituting the form (B ll) into the input/output relation (B8) gives the complex 
modulation u s (£) of the returned pulse. When this is substituted into (B6) one obtains: 


E s (<) = x|,4|/2 {/(*'_) cos [<£(£'_) - C Jot- 4- 0.4 - 'i’K 
/(£+) cos 4* 0.4 4- ^7T 


4- 

}> 


(B 13) 


/ jftr 

t- =t 

C 


2 a 

4 > 

c 


(B 14) 


where 

2a t' — t ^ r 

c ^ c 

a is the radius of the oval, and <p A is the phase angle of .4. Equations (B 13) and (B 14) 
have the following interpretation. The returned pulse can be thought of as consisting of 
the incident pulse reflected twice, each time with reflection coefficient .4/2 and phase 
shifts 4> a = ±37r/4. The two reflections occur at z — ±a. The time shift -R r /c accounts 
for the travel time from the point z = 0 to the receiver while the additional time shifts 
±2a/c account for the round-trip time required to travel from z - 0 and the two points 
z — ±a. 

For the oval with the slice of entrained ambient fluid, the behavior of H(uj) was more 
complicated than (B 11). In this case the procedure followed was to first calculate H(oj) 
for us spanning an interval Acj consisting of few pulse bandwidths about the center 
frequency. Actually a convenient auxiliary quantity T(o/) was calculated such that 

T( u)e tkRr 

V ' (B 15) 




Rr 


Next T(cj) was multiplied by a Hanning window and Fourier analyzed such that: 


T ^) - ^2 Tmex p 


2?r ( \ 
i — — miijj — u Jn) 

Ac o 


(B 16) 


Knowing the coefficients T m one can obtain the complex modulation of the received signal 
by substituting (B 16) and (B 15) into (B8): 




(B 17) 


Rectangular and Gaussian forms of unit amplitude were implemented for the amplitude 
modulation: 


m 


={°: 


m 




(B 18) 


1*1 5% r/2, 
elsewhere ’ 

where r is the pulse duration. In this paper only the results for the Gaussian pulse are 
plotted since the rectangular pulse gave very close RCS results (within 0.25 db for a few 
spot checks). For FM a linear frequency variation (chirp) was considered: 

<p{t) = 7vat 2 . (B 19) 

Within the pulse width, the spread A/fm in the local frequency is 

A/fm = |a|r. (B20) 

The sign of a (up-chirp vs. down-chirp) was found to not matter. Table 5.2.1 gives r and 
A/fm for Gilson’s tests. 
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The instantaneous power flux density vector associated with the incident pulse (B3) 
and (B 4) is 

Si = E; x Hi = zc£ 0 f 2 (t)cos 2 (<t>(t) -u 0 1) . (B21) 

Finally, the instantaneous scattered power flux density vector associated with (B6) 
and (B 7) is 

S s (t) = zce 0 |E 9 (t)| 2 . (B 22) 

From these we calculate an instantaneous RCS. 

RCS(t) = (B23) 
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